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Abstract 

We consider an Yb-doped double-clad fiber laser in a unidirectional ring 
cavity containing a polarizer placed between two half- wave plates. Depending 
on the orientation of the phase plates, the laser operates in continuous, Q- 
switch, mode-lock or unstable self-pulsing regime. An experimental study 
of the stability of the mode locking regime is realized versus the orientation 
of the half-wave plates. A model for the stability of self-mode-locking and 
cw operation is developed starting from two coupled nonlinear Schrodinger 
equations in a gain medium. The model is reduced to a master equation in 
which the coefficients are explicitly dependent on the orientation angles of 
the phase plates. Analytical solutions are given together with their stability 
versus the angles. 



1 Introduction 



Self-mode-locking has been successfully obtained in both low power and 
double-clad fiber lasers using Additive-Pulse Mode locking technique through 
nonlinear polarization rotation PQ-[E]- The experimental configuration gen- 
erally consists in a ring cavity containing a polarizing isolator placed between 
two polarization controllers. The basic principle is the following. The polar- 
ization state evolves nonlinearly during the propagation of the pulse in the 
fiber due to the combined effects of self-phase modulation and cross-phase 
modulation induced on the two orthogonal polarization components, both 
resulting from the optical Kerr effect. A polarization controller is adjusted 
at the output of the fiber such that the polarizing isolator lets pass the cen- 
tral intense part of the pulse but blocks the low-intensity pulse wings. This 
technique of nonlinear polarization rotation for passive mode locking leads 
to stable, self-starting pulse trains. Theoretical approaches to this prob- 
lem can be divided into two categories. The first one has been proposed 
by Haus et al j7] and consists in writing a phenomenological scalar equation 
(the master mode-locking equation) assuming that all effects per pass are 
small. The model includes the group velocity dispersion (GVD), the optical 
Kerr nonlinearity and a gain medium. No birefringence is taken into ac- 
count. Although the mode-locking properties of the laser can be described 
through this approach for positive and negative GVD, the stability of the 
mode-locked solutions as a function of the orientation of the polarization 
controllers cannot be described in this way. On the other hand Kim et al 
[S] use an approach based on two coupled nonlinear Schrodinger equations. 
The medium is assumed to be birefringent. The GVD together with the Kerr 
effect are taken into account. The orientation angles of the eigenaxis of the 
fiber at both sides of the polarizer are explicitly taken into account through 
a periodic perturbation. Numerical simulations show that mode-locking can 
be achieved for negative GVD. To our best knowledge, none of the models 
take into account all the characteristics of the fiber (birefringence, GVD, 
gain, optical Kerr effect) and the action of the orientation of the polarization 
controllers. 

The aim of this paper is to investigate experimentally and theoretically 
the mode-locking properties of the Yb-doped double-clad fiber laser. Mode- 
locking is achieved through nonlinear polarization rotation in a unidirectional 
ring cavity containing a polarizer placed between two half-wave plates. This 
paper is organized as follows. Section II is devoted to the experimental re- 
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suits. We first briefly recall the operating regime of the laser as a function 
of the orientation of the two half-wave plates. A stability diagram of the 
mode-locked laser is presented fSj. The parameters required in the mod- 
elling section are measured. The birefringence of the double-clad fiber is 
first determined using the magneto-optic method ^U]. The total dispersion 
has been estimated from other experiments where compression of the pulses 
was realized I n section |3] the stability of the mode-locking regime 

is theoretically investigated. The analysis is based on the coupling between 
the two linearly polarized eigenstates of the fiber jHUSl- The propagation 
problem is modelled through two nonlinear Schrodinger equations in a gain 
medium where the GVD, the nonlinear coefficient and the gain are assumed 
small over the cavity length. The eigenaxis of the fiber at each side of the 
polarizer are referenced with respect to the polarizing axis. These angles can 
be adjusted with two half-wave plates and will thus be considered as variable 
in the analysis. The polarizer is treated as a periodic perturbation acting 
as a projecting operator. Using the expressions derived in the propagation 
problem, we calculate the electric field amplitude just after the polarizer as 
a function of the corresponding amplitude at the previous round-trip. The 
round-trip number is then considered as a continuous variable, which allows 
us to use a multiscale analysis approach. The final equation is a complex 
Ginzburg-Landau equation whose coefficients explicitly involve the orienta- 
tion angles of the half-wave plates. Explicit solutions are given and their 
stabilities are discussed. A stability diagram of the mode-locked regime is 
calculated and compared to the experimental diagram. 



2 Experimental results 

In this section we briefly recall the results obtained with an Yb-doped double- 
clad fiber laser in a unidirectional ring cavity operating at 1.08 fim 0. A 
polarizing isolator is placed between two half-wave plates. Although this con- 
figuration is known to generate ultrashort pulses in low power fiber lasers, 
it can also lead to very different regimes depending on the orientation of 
the polarization elements, such as cw operation, Q-switching, mode-locking, 
or unstable emission regime. We thus perform a systematic study of the 
operating regime as a function of the relative orientation of the two half- 
wave plates. The experimental setup is shown in figure H The 4 meter long 
Yb-doped double-clad fiber is side-pumped using the v-groove technique j3] . 



2 



Yb-doped 
double-clad fiber 



V-groove 




Pump 975 nm 



A/2 



A/2 



output 



Optical 
Isolator 



Figure 1: Experimental setup. 



The launching efficiency is about 70 percent. This pumping scheme is very 
convenient for ring laser applications because the two fiber ends remain free. 
The fiber core diameter is 7 /zm allowing a single-mode laser signal propaga- 
tion. The inner cladding is a 125 x 125 pm 2 square leading to a multimode 
pump propagation. The geometry of the inner cladding ensures an efficient 
pump absorption. Indeed, the pump power is almost completely absorbed 
along the doped fiber. In addition, multimode propagation is required for 
high injection efficiency of large area laser diodes. Two single-mode fibers 
at 1 /iin are spliced at both ends of the double-clad fiber leading to a total 
fiber length of about 8.5 m, the cavity length is about 9 m. Fiber ends are 
angle-cleaved in order to avoid any Fresnel reflection which could generate 
a standing- wave operating regime. The splicing losses are less than 0.1 dB. 
The doped fiber is pumped at 975 nm with a 4 W semiconductor laser. A 
bulk polarizing optical isolator is used to obtain a travelling wave laser. The 
latter is placed between two half- wave plates at 1.08 /im. Each plate can be 
rotated in a plane perpendicular to the propagation direction. A partially 
reflecting mirror ensures an output coupling of about 90 percent. The output 
beam is analyzed through different means. A fast photodiode is used when 
the laser is cw, Q-switched or unstable. In the regular mode-locking regime, 
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we use an optical autocorrelator to measure the duration of the pulses. Note 
that there is no compression system in the cavity. An optical spectrum ana- 
lyzer (Advantest) together with a microwave spectrum analyzer (Tektronix) 
are also used when needed. Finally, a powermeter allows us to measure the 
average output power. 

We investigate now the operating regime of the laser as a function of the 
orientation of the two waveplates. 61 (6 2 ) is the orientation of one eigenaxis 
of the first (second) half-wave plate referenced to the input (output) polarizer 
of the optical isolator. The angle between the input and output polarizer of 
the isolator is 45 degrees. For a given pumping power we record the evolution 
of the output intensity versus 6\ and 62- Practically, we fix 6\ and measure 
the particular values of 62 for which a change in the operating regime occurs. 
It is convenient to represent these results in the plane (61,62)- It allows 
to easily identify the output regime of the laser. The experimental results, 
obtained for a pumping power of about 2.75 W, are summarized in figure 
121 We have experimentally observed that the different output regimes are 




30 60 90 120 150 180 

9, (deg) 



Figure 2: Operating regime of the laser as a function of the orientation of the two half- 
wave plates 9i and 02 ■ The pump power was 2.75 W. ML stands for mode- locked regime 
and cw for continuous wave regime. 

periodic versus 6i and 6 2 with a period of 90 degrees. Figure El shows that 
the laser essentially operates in cw, mode-lock or unstable regime while the 
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Q-switch behavior appears near the boundary between cw and mode-lock 
regimes. Note that since the operating regime depends on the pumping rate, 
the mapping of figure 121 is pump power dependent. Typically, when the pump 
power decreases (increases), the ranges of mode-locking decrease (increase). 

The different dynamical behaviors observed in our experiment have been 
discussed in details in reference 0. Here, we concentrate our attention on 
the mode- locked operating regime. In figure E3 we give the results obtained for 
Q\ = 20 deg and 9 2 = deg. In this case, the laser is first cw just above thresh- 
old (figure Ek). While the pump power is increased, an irregular self-pulsing 
regime appears which finally becomes a regular self-mode-locking regime. As 
previously discussed in the introduction, the basic principle responsible for 
the mode-locking is well known in such configurations. It is based on Kerr 
effect through nonlinear polarization rotation. This technique is also called 
Polarization- Additive Pulse Mode- locking. The repetition rate of the pulses 
corresponds to the free spectral range of the cavity, about 20 MHz in our 
case. Auto-correlation traces showed a coherence spike of about 150 fs on a 
longer pedestal of about 60 ps jBj. A typical spectrum in the mode-locking 
range is given in figure The mode-locked spectrum bandwidth is more 
than 25 nm. This spectrum suggests that shorter pulses can be achieved by 
a suitable dispersion compensation. We have recently shown that femtosec- 
ond pulses are obtained with a grating pair inserted in the cavity jSj. In the 
present mode-locked regime, we measure 600 mW output power for 3.25 W 
pump power, which corresponds to high energy pulses of about 30 nJ. 

Consider now the measurements of parameters required in the theoretical 
section. We first measure the beating length Lb of the double-clad fiber using 
the magneto-optic method JUj- The method consists to apply a sinusoidal 
magnetic field produced by a coil in a small portion of the fiber and to move 
the coil along the fiber axis. The fiber is placed between an analyzer and a 
polarizer. This system is then illuminated and analyzed in transmission. The 
detected output power is a periodic signal versus the coil displacement. The 
period in the spatial domain coincides with the beating length of the fiber. 
We found Lb = 50 cm leading to a birefringence parameter K = l/(Lg) = 
2 m _1 . We must also consider the birefringence of the undoped fibers. It 
has been measured to be of about 1 m . In the following we will consider 
an average value of the birefringence of about 1.5 m _1 . A second important 
parameter is the GVD. It has been estimated from compression experiments 
[TJ]. We find (3 2 = 0.026 ps 2 /m. 
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Figure 3: 6i = 20 deg and 62 — Odeg. (a) laser characteristic and (b) typical optical 
spectrum in the mode-locked regime. 



3 Theory 

The aim of this section is to develop a theoretical model able to describe the 
mode-locking properties of the Yb-doped fiber laser as a function of the ori- 
entation of the two half-wave plates. For that, we consider a birefringent fiber 
in a unidirectional ring cavity with an intracavity polarizer. The orientation 
of the eigenaxis of the fiber at each sides of the polarizer are unspecified and 
will be taken as variable parameters. This is equivalent to considering that 
the eigenaxis of the fiber are fixed and two half- wave plates are used to adjust 
their orientations with respect to the polarizer. 
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3.1 Derivation of a master equation 
3.1.1 Propagation in the fiber 

The starting point are the equations giving the evolution of the two polar- 
ization components in a gain medium with Kerr nonlinearity and GVD. In 
the framework of the eigenaxis of the birefringent fiber moving at the group 
velocity, the pulse envelope evolution is described by the following system 

EM 

id z u -Ku- y<9 2 w + 7 (u \u\ 2 + Au \v\ 2 + Bv 2 u*) = ig (l + -^d 2 J u, (1) 

id z v + Kv- ^d 2 v + 7 (v \v\ 2 + Av \u\ 2 + Bu 2 v*) = ig ^1 + -^d$J v, (2) 

where d t denotes the partial derivative operator J^. /3 2 (in ps 2 /m) is the 
GVD coefficient. K (in m _1 ) is the birefringent parameter and is related 
to the x and y refractive indexes through the relation K = n(n x — n y )/X, 
where lambda is the optical wavelength. 7 = 27rn 2 /(AA e //) is the nonlinear 
coefficient, n 2 (in m 2 /W) is the nonlinear index coefficient and A e ff (in m 2 ) 
the effective core area of the fiber. A and B are the dielectric coefficients. 
In isotropic media, A = 2/3 and 5 = 1/3 [T^j. This is the case with our 
silicate fibers, g is the gain coefficient (in m _1 ) and u g is the spectral gain 
bandwidth (in ps _1 ). At first order the gain coefficient is fixed by the fact 
that it compensates the losses. Note that polarization mode dispersion is 
not taken into account since the cavity length is short (about 9 m). For 
numerical simulations we will use the following values for the parameters: 
K = 1.5 m" 1 , (3 2 = 0.026 ps 2 /m, L = 9 m, 7 = 3 ■ 10~ 3 W^m" 1 and 
Ug = 10 13 s -1 . 

We now assume that the GVD /3 2 , the nonlinear coefficient 7, the gain 
filtering p = g/u 2 are small over one round-trip of the cavity. We introduce 
a small parameter e and replace these quantities by e/3 2 , £7 and ep. We then 
look for solutions of the system (JTJ2J) under the form of a power series in e, 

as 

u = u + eui + O (e 2 ) , (3) 
v = v + evi + O (e 2 ) . (4) 

At order e° it yields 

u = u e^ iK >, (5) 
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v = v e 



(g+iK)z 



(6) 



where u and v do not depend on z, ie depend on the time variable t only. 
Making the transformation 



Ul = u 1 (z,t)e (9 ~ lK)z , 
v 1 = d 1 (z,t)e^ + ^ z , 



(7) 
(8) 



the equations obtained at order e can be integrated with regard to z to yield, 
for the component u: 



e 2gz _ y 

+ry (u \u \ 2 + Au \v \ 2 ) 



e (2g+iiK)z _ 1 ( 9 ) 

2g ' "'~" u "° 2g + UK 
To get the complete expression for u(z), we notice that 

u(0) = u + e«i(0) + O (e 2 ) . 
Then the two components of the wave amplitude can be written as 



(10) 



u(z) = u(0)e (9 ~ iK)z + e z ( p - ^ ) d 2 u{Q) 



+17 (u(0) |«(0)| 2 + Au(0) K0)r) 



2\ e 



2gz _ y 



+i7Bv(Q) 2 u(Q)* 
v(z)= v(0)e {9+tK)z + e 



2g 

e (2g+AiK)z _ i 



:n) 



z\p- 



2g + UK 



+i7 (v(0) \v(0)\ 2 + Av(0) \u(0)\ 2 ) 



d{v(0) 

e 2gz _ i 



+i-fBu(0) 2 v(0y 



2g 

e {2g-4iK)z _ l 

2g - UK 



(12) 



,GH-uO* + o (e 2 ) . 



3.1.2 Polarizer 



Figure 0] shows the respective orientation of the eigenaxis of the fiber before 
and after the polarizer (optical isolator). (u-,V-) are the field components 
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just before the polarizer, (u + , v + ) just after it. The effect of the polarizer can 
be written as 

£M sm9 - ) ("-")■ (13) 

where (3 is its transmission coefficient ((3 < 1). Immediately after the po- 
larizer, the field has a well-defined linear polarization. Denote by f n its 
amplitude at the beginning of the n th round trip. Equation (JT2|1 can be 
written as 

/ n+1 = /3(cos0_ sin0_)(^-' n Y (14) 



v. 



and 

U+,n+l \ ( COS0+ 



I - I SHI a I fn+U ( 15 ) 

where (u n ,v n ) are the field components during the n th round. The field 
components u +n = u n (0) and v +n = v n (0), at the entrance of the fiber 
immediately after the polarizer, are transformed respectively into u n (L) = 
u_ n and v n (L) = v_ n , at the end of the fiber immediately before the polarizer 
(L is the fiber length). 

w n (0) and v n (0) are proportional to f n . Using the propagation formulas 
(fTTjl and (fP2]l. u n (L) and v n (L) can be computed. Then / n+1 is computed as 
a function of f n using (JHJ). This yields 

/„+l=/3e» £ {<?/„+£ (p-^^)LQ^/n + ^/n|/n| 2 }+0(£ 2 ), (16) 
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with 
and 



P = 7 



Q = cos 

e 2gL _ j 



6L) cos - i cos(9 + + 6L) sin KL, 



(17) 



2<? 



cos(# + — 6L) cos i^L — i cos( 



6J)smKL 



A 



■ sin 29+ [sin(fl+ + 0_) cosKL -isin(9+ - 0_) sini^L] 



+ — sin 2# H 



sin 0. cos9.e~ iKL 



+ cos 0+ sin 6Le 



e (2g+4iX)L _ 1 
e (2g-4iK)£ _ 1 



(18) 



2g - 4iiT 

3.1.3 Gain threshold and continuous limit for f n 

The dominant behavior of (f n ) is given by the zero order term in (|16j): 

/ n+1 = /3e 9L Q/ n + O (e) . (19) 

A steady state over a large number of round trips could be reached if this 
first order evolution is the multiplication by a phase factor, but only in this 
case. Therefore the modulus of f3e 9L Q must be 1. This gives the threshold 
gain value: 



9o = ^ln(/3 2 |Q| 2 ) 



2L 



In ((3 2 [cos 2 (9+ - 9 J) - sin 29 + sin 29 _ sin 2 KL] ) . 



(20) 



Denote by e 1 ^ the quantity f3e 9oL Q. The exact condition is |/5e 9L Q| = 1 + 
O (e), the gain g therefore writes g = go + £9i, hi which g>i is still free. 
Expanding e £9lL in a power series of e, equation (fTo^) becomes 



= e** (1 + e^L) /„ + e(^p - -f J Le^/ B + i£ -Q~fn \fn\ 2 + O (e 2 ) . 

(21) 

Consider now some function / of a continuous variable 2, obeying an 
equation of the form 



idj = {A + ieB) f + eCd 2 f + eVf \ff 



(22) 
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where A and B are real, and C and T> complex coefficients. We look for 
solutions under the form / = f + efi + . . .. The solution at first order is 



/o = f e- iAz . (23) 

We write 

h=~h{z,t)e- iM , (24) 

and obtain: 

f(L) = [/(0) + e (B/(0) - ^ t 2 /(0) - fl>/(0) l/(0)| 2 ) e^ L + O (^ 2 ) . 

(25) 

Equations (J2T]) and (|25|) are identified according to /(0) = / n and /(£) = 
/n+ijto yield 

e -iAL = e # ) (26) 
B = 0i, (27) 
C = f + ip, (28) 



Then the continuous equation (j22|) yields some interpolation of the dis- 
crete sequence (f n )- The continuous approximation is relevant when the 
number of round trips is very large; further, mode-locked pulses actually cor- 
respond to this situation. The small correction of order e in equations (|22j) 
or (|25|) gives account for the variations of / on 'propagation distances', which 
are here number of round trips, very large, of order 1/e. This can be shown 
in a rigorous way using the multiscale formalism, commonly used for the 
derivation of the model equations in the soliton theory ^Hj- We introduce a 
slow variable ( = ez, in such a way that 

d z = d £ + ed c . (30) 

The values of ( about 1 correspond to number of round trips about 1/e. f 
is expanded as above as fo + efi + . . ., and the first order obviously yields 
we make the transform (J2*lj) and f\ must satisfy the following equation: 



^^ c fo + id & h = iBfo + C%f + Vf f * . (31) 
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Equation (|31|) can be written as 

dJi=Hfo). (32) 

The long distance evolution of fo is obtained from the requirement that j\ 
does not grow linearly with z. This yields d$fi — and 

idjo = iB~h + Cd 2 J + Vf |/ | 2 • (33) 
Or, using real coefficients only: 

idjo = igJo + + ip) d tfo + (V + Wi) f |/o| 2 • (34) 

Equation or is the cubic complex Ginzburg-Landau (CGL) equa- 
tion. Note that equation (|3*3*|l is formally identical to the master equation 
proposed by Haus et a/ jjj . However, the coefficients of (|3*3*|l explicitly depend 
on the orientation of the eigenaxis of the fiber at both sides of the polarizer. 
An essential feature is the arising of a nonlinear gain or absorption Z>j, which 
results from the combined effects of the nonlinear rotation of the polariza- 
tion, the losses due to the polarizer, and the linear gain. The value and the 
sign of T>i depend on the angles 9 + and 6L between the eigenaxis of the fiber 
and the polarizer. 

3.2 Different regimes for the CGL equation 

3.2.1 The stationary solution and its modulational instability 

The CGL equation (J34|) admits the following nonzero stationary solutions, ie 
with a constant modulus: 

/i = Ae i{< - nt) (35) 

with 

il 2 = - (V t \A\ 2 + 9l ) (36) 
P 

and 

K = ^(V l \A\ 2 + g 1 )-V r \A\ 2 . (37) 

In the particular case Q = it is a constant solution (independent of t), with 
a fixed amplitude 

A = M, (38) 
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and 

« = ^91- (39) 



It exists only if T>ig\ < 0, ie if the excess of linear gain g\ and the effective 
nonlinear gain T>i have opposite signs. We perform a linear stability analysis 
for the constant solution f c = Ae iK< *. We seek for solutions of the form 
fi — fc (1 + u )i u being very small. It is seen that it must satisfy the following 
equation: 

-ku + id ( u = ig x u + + ipj d t u + VA ^ O + u *) • ( 40 ) 

Due to the term u*, equation (j4*Uj) is not linear. The eventual instability is of 
the type discovered first by Benjamin and Feir in the frame of water waves 
[To] . It has been shown that this kind of instability involves the nonlinear 
interaction of two modulation terms with conjugated phases with the square 
A 2 of the constant solution [T2j. Thus we write u as: 

u = Ul e {K ~ iujt) + u 2 e {x ^ +iuJt) (41) 

It is found that A satisfies an equation of the form 

A 2 + 2b\ + c = 0. (42) 



Using the expression (}38|) of the fixed amplitude A, we get the following 
values of the real constants b and c: 

b = 9l + pco 2 , (43) 



^+H 2 +* 2 -4 + !>- (44 » 



A short analysis of equation ()42j) show that a necessary condition for u to re- 
main bounded is that b is positive. This proves that modulational instability 
occurs as soon as the excess of linear gain g\ is negative and the nonlinear 
gain T>i is positive. Recall indeed that the existence of the constant nonzero 
solution requires that these two quantities have opposite signs. When g x > 
and T>i < 0, we show that the modulational instability never occurs. Numer- 
ical computation shows indeed that the quantity (P r /3 2 /2 +Vip) is always 
negative, which ensures, using elementary analysis, that c is always positive, 
and the result follows. 
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If we admit that the excess of linear gain g\ will self-adjust to a value for 
which a stable solution exists, the above discussion shows that the stability 
of the constant solution depends on the sign of the effective nonlinear gain 
T>i. It is stable when T>i < only. We expect to observe continuous laser 
emission when the constant solution is stable, and only in this case. Thus 
the domains of continuous emission should coincide with the regions where 
T>i is negative. The sign of T>i as function of the angles (^r#_, \6 + ) is drawn 
on figure El It is easily checked that the relation between the angles 9\, 6 2 
used in the experiments and the angles 9 + , 6L defined by figure |U is such 
that 6>i = -^0- + 6*io and 6 2 = \0 + + 9 2 o, with some fixed value of 6\q and 
#20- The white domain on figure corresponds to the negative values of 




=^6_ (deg) 

Figure 5: Sign of the effective nonlinear gain T>i as a function of the orientation of the 
polarizer. White: X\ < 0, and the constant solution of the CGL equation is stable. Gray: 
T>i > 0, and the constant solution of the CGL equation is unstable. 

the nonlinear gain T>i. It is thus the domain where the constant solution of 
the CGL equation is stable. Continuous laser emission should occur in this 
domain. 
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3.2.2 Localized solutions 



A localized analytical solution of the CGL equation (J3*lj) can also be written 
[THj . It has the following expression: 

f ± = a(t)e i{dXna{t) -^\ (45) 

with 

a(t) = BCsech (Bt) (46) 

and 



B - y pd* - 3 p - w (47) 



C= . /; *M + ® (48, 



2 (j3 2 Vi - 2pV r ) ' 



3 (p 2 V r + 2pVi) + \/9 (f3 2 V r + 2pVif + 8 {faVi - 2pV r f 
2((3 2 V i -2pV r ) ' 



(49) 



- 9l (4pd + (3 2 d 2 - (3 2 ) 
U ~ 2(pd*-p-f3 2 d) ■ (50) 

The inverse B of the pulse length is real only if the quantity T = (pd 2 — p — (3 2 d) 
and the excess of linear gain g\ have the same sign. The regions where T 
is either positive or negative are specified on figure El The background with 
zero amplitude is stable when the excess of linear gain gi is negative, and 
unstable in the opposite case. When gi > 0, the exact localized solution 
()45|) is unstable due to the instability of the background [T%] . Qualitatively, 
localized pulse formation can be expected when the excess of linear gain gi 
is negative and the nonlinear gain T>i is positive, as is suggested on figure [7| 
The effective self-phase modulation T> T is always negative, and the dispersion 
/?2 is positive. Therefore their conjugated effect leads to the increase of the 
pulse width. If g\ is positive and T>i negative (figure [7^,), the nonconservative 
effects decrease the amplitude at the top of the pulse, and increase it at the 
bottom: no stable localized pulse can be formed. If on the contrary g\ is 
negative and T>i > (figure |7b), energy appears at the top of the pulse and 
disappears at its bottom, yielding some pulse narrowing, which could be ex- 
pected to balance the broadening caused by the nonlinear index variations. 
In fact the analytical localized solution (|43|) is never stable when j3 2 V r < 0, 
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Figure 6: Sign of the excess of nonlinear gain value g\ for which the exact localized 
solution exists, as a function of the orientation of the polarizer. White: T and g\ positive. 
Gray: T and g\ negative. 




Figure 7: Schematic representation of the effect of the nonlinear gain X>j, of the excess 
of linear gain g\ , and of both the dispersion /?2 and the effective self-phase modulation £>,. 
on a localized pulse. 



but can be stabilized when higher order nonlinear terms are taken into ac- 
count, as shows the study of the so-called quintic CGL equation [THj. The 
stabilizing term should be a quintic nonlinear absorption. Observe that here 
the effective nonlinear gain T>i follows from the linear gain and the nonlinear 
evolution of the polarization in the fiber, without considering cubic nonlinear 
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gain in the initial model (0I2J)- In an analogous way, an effective quintic non- 
linear gain term in the master equation (|33|) does not require that a quintic 
nonlinear gain term is present in the initial model, either a cubic nonlinear 
gain or a quintic nonlinear index should be enough. (JTJ2J). The precise deter- 
mination of this higher order term is left for further study. On these grounds, 
the exact solution (J45j) corresponds to a potentially stable pulse when Q\ < 
and T>i > 0, and to a completely unstable pulse when gi > and T>i < 0. 

3.2.3 The different regimes 

According to the conclusions of the previous section, and assuming that 
the excess of linear gain g\ will self-adjust to the value for which a stable 
solution exists, mode-locked laser emission can be expected in the region 
where Pj > and T < 0. This is the domain in dark gray marked ML 
on figure |S1 In the region where Z>, is negative, the constant solution of 
the CGL equation is stable, and continuous laser emission is expected to 
occur. It is the white domain marked cw on the figure |H1 In the region where 
T is positive and T>% negative, no stable localized neither constant solution 
exist. The laser behavior is expected to be unstable in this region, in light 
gray on the figure |H1 The theoretical results summarized on the figure |H1 
can be compared to the experimental results of figure EJ A discrepancy 
between the two figures appears at first glance: the periodicity of the behavior 
regarding the variable 9<i or 9 + /2 seems to be 45 degrees according to the 
theory, while it is 90 degrees according to the experiment. Indeed, the two 
elongated domains of mode-locked behavior are not equivalent according to 
experimental observation, while they seem to be identical on the theoretical 
figure |H1 A contour plot of the effective nonlinear gain T>i as a function of 
the angles is drawn on figure EJ It is seen that the values of T>i are not the 
same in the two theoretical domains of mode-locked emission. Thus they 
are in fact not equivalent, and the periodicity of the theoretical result is 90 
degrees, as that of the experimental one. Further it has been shown [TH] that 
an excessive value of the nonlinear gain might prevent pulse stabilization. A 
more accurate analysis would thus very likely show that the regions where 
T>i takes its highest values are out of the domain of stability of the localized 
pulse. This confirms that the two horizontal domains of mode-locked regime 
drawn on figure |S1 are not equivalent. Further, it could explain a part of the 
discrepancy between theoretical and experimental results. 

Before we conclude, let us consider the influence of the group velocity 
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dispersion on the stability of the mode-locking solutions. This is important 
because GVD compensation is required to generate subpicosecond pulses. 
Experimentally, we have obtain 666 fs pulses with a grating pair inserted 
in the cavity jHHHJ. Under these conditions, the total GVD is about /3 2 = 
0.0005 ps 2 /m. Hence, stable ultrashort pulse generation is possible for such 
low values of the GVD. We have tested our master equation for decreasing 
values of /?2- We have found that a lowering of P2 results in a reduction 
of the stable mode-locking regions. In addition, mode-locking completely 
disappears for 02 < 0.0015 ps 2 /m. These results are in agreement with the 
predictions of the Haus's model [7j but not with our experimental results. 
The validity of the model presented here is therefore limited to GVD values 
not too close to zero. Higher order terms of the GVD could be required to 
improve the model near the zero-dispersion point. This problem deserves 
further experimental and theoretical work. 
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Figure 9: Contour plot of the effective nonlinear gain X\ as a function of the angles 0+ 
and 9- between the polarizer and the eigenaxis of the fiber. 



4 Conclusion 

We have given a rigorous derivation of the master equation, a CGL equa- 
tion, that allows to discuss theoretically the different operating regimes of 
the passively mode-locked Yb-doped double-clad fiber laser. This equation 
involves a single amplitude, but gives an account of the evolution of the po- 
larization inside the fiber through its coefficients. An effective nonlinear gain 
arises from the conjugated effects of the linear gain and the self-phase mod- 
ulation. Its sign depends on the orientation of the polarizer. The existence 
and stability of the constant nonzero solution of the CGL equation has been 
discussed. It gives an account of the continuous emission regime of the laser. 
An analytical localized solution has also been given. Although it is properly 
speaking unstable, it allows to discuss the formation of localized pulses, ie the 
mode-locked regime of the laser. An unstable regime has also been identified. 
These theoretical results are in good accordance with the experiments. The 
study of higher order corrections to this model should improve these results. 



19 



References 



[1] K. Tamura, H.A. Haus and E.P. Ippen, Self-starting additive 
pulse mode-locking erbium fiber ring laser, Electron. Lett. 28, 
2226-2228, 1992. 

[2] H.A. Haus, K. Tamura, L.E. Nelson and E.P. Ippen, Stretched- 
pulse additive pulse mode-locking in fiber ring lasers: theory and 
experiment, IEEE Jour. Quant. Electron. 31, 591-598, 1995. 

[3] V.J. Matsas, T.P. Newson, D.J. Richardson and D.N. Payne, 
Selfstarting passively mode-locked finer ring soliton laser ex- 
ploiting nonlinear polarization rotation, Electron. Lett. 28, 
1391-1393, 1992. 

[4] V. Cautaerts, D.J. Richardson, R. Pashotta and D.C. Hanna, 
Stretched pulse Yb3+:silica fiber laser, Opt. Lett. 22, 316-318, 
1997. 

[5] R. Hofer, M. Hofer and G.A. Reider, High energy, sub- 
picosecond pulses from a Nd-doped double-clad fiber laser, Op- 
tics Com. 169, 135-139, 1999. 

[6] A. Hideur, T. Chartier, M. Brunei, S. Louis, C. Ozkul and F. 
Sanchez, Generation of high energy femtosecond pulses from 
a side-pumped cYb-doped double-clad fiber laser, Appl. Phys. 
Lett. 79, 3389 (2001). 

[7] H.A. Haus, J.G. Fujimoto and E.P. Ippen, Structures for ad- 
ditive pulse mode locking, J. Opt. Soc. Am. B 8, 2068-2076, 
1991. 

[8] A.D. Kim, J.N. Kutz and D.J. Muraki, Pulse-train uniformity 
in optical fiber lasers passively mode-locked by nonlinear po- 
larization rotation, IEEE Jour. Quant. Electron. 36, 465-471, 
2000. 

[9] A.Hideur, T.Chartier, M.Brunei, M.Salhi, C.Ozkul and 
F.Sanchez, Mode-lock, Q-switch and cw operation of an Yb- 
doped double-clad fiber ring laser, Optics Com. 198, 141-146, 
2001. 



20 



[10] T. Chartier, A. Hideur, C. Ozkul, F. Sanchez, G. Stephan, 
Polarization properties of single-mode optical fibers with ran- 
dom fluctuations of birefringence, Appl. Optics. 40, 5343-5353 
(2001). 

[11] M.L. Dennis and I.N. Duling III, Experimental study of side- 
band generation in femtosecond fiber lasers, IEEE Jour. Quant. 
Electron. 30, 1469-1477, 1994. 

[12] A. Hideur, T. Chartier, C. Ozkul, M. Brunei and F. Sanchez, 
Experimental study of pulse compression in a side-pumped Yb- 
doped double-clad mode-locked fiber laser, Appl. Phys. B 74, 
121-124 (2002). 

[13] G.P. Agrawal, Nonlinear Fiber Optics, Academic Press, Second 
Edition, 1995. 

[14] D.J. Ripin and L. Goldberg, High efficiency side-coupling of 
light into optical fibers using embedded V-grooves, Electron. 
Lett. 31, 2204-2205, 1995. 

[15] T. Taniuti and C.-C. Wei, Reductive perturbation method in 
nonlinear wave propagation I. J. Phys. Soc. Japan, 24, 941-946 
(1968). T. Taniuti and N. Yajima, Perturbation method for a 
nonlinear wave modulation. Ill J. Math. Phys., 14 (10), 1389- 
1397 (1973). 

[16] T.B. Benjamin and J.E. Feir, The desintegration of wave trains 
in deep water. Part.l. Theory, J. Fluid. Mech. 27 (3), 417-430 
(1967). 

[17] J.T. Stuart, F.R.S. and R.C. DiPrima, The Eckhaus and 
Benjamin-Feir resonance mechanisms, Proc. R. Soc. Lond. A 
362, 27-41 (1978). 

[18] J.M. Soto-Crespo, N.N. Akhmediev, V.V. Afanasjev, and 
S.Wabnitz, Pulse solutions of the cubic-quintic complex 
Ginzburg-Landau equation in the case of normal dispersion, 
Phys. Rev. E 55 (4), 4783-4796 (1997). 



21 



